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We formulate a new method to calculate the gravitational reaction force on a particle of mass 
/i orbiting a massive black hole of mass M. In this formalism, the tail part of the retarded Green 
function, which is responsible for the reaction force, is calculated at the level of the Teukolsky 
equation. Our method naturally allows a systematic post-Minkowskian (PM) expansion of the tail 
part at short distances. As a first step, we consider the case of a Schwarzschild black hole and 
explicitly calculate the first post-Newtonian (1PN) tail part of the Green function. There are, 
however, a couple of issues to be resolved before explicitly evaluating the reaction force by applying 
the present method. We discuss possible resolutions of these issues. 



I. INTRODUCTION 



In the black hole perturbation approach to gravitational radiation from an inspiralling compact binary, one calculates 
the gravitational waves emitted by a point particle orbiting a black hole, assuming the mass of the particle (fi) is 
much less than the mass of the black hole (M). At the lowest order in n/M, the particle moves along a geodesic 
trajectory in the black hole geometry. Already at this lowest order, this approach has been proved to be very powerful 
for evaluating general relativistic corrections to the gravitational wave forms from a compact binary, even for neutron 
star- neutron star (NS-NS) binaries jij. However, in such calculations one inevitably has to assume that the radiation 
reaction time scale is sufficiently long compared to the orbital period, so that the adiabatic approximation to the 
orbital evolution is valid. That is, the orbital evolution is determined by the energy and angular momentum balance 
between the orbit and the graviational radiation. 

However, there are several situations in which such an approximation cannot be justified. For example, when the 
particle is nearly in the last stable circular orbit, the radiation reaction time scale becomes comparable to the orbital 
period. When the particle is in a very eccentric orbit or in a scattering orbit, presumably one has to determine the 
orbital evolution not by averaging over many periods but by each scattering event. 

Another case is when the massive object is in fact a rotating black hole. In this case the geometry is Kerr, and the 
geodesic orbit is determined not only by the energy and the z-component of the angular momentum but also by the 
Carter constant [||. Since the Carter constant is not associated with the Killing vector field of the Kerr geometry, 
its evolution cannot be determined by evaluating the gravitational radiation at infinity, even though the adiabatic 
approximation may still be valid. 

In order to treat these situations, it is therefore necessary to evaluate the correction to the equation of motion of a 
particle in the background geometry. That is, we must know the reaction force of 0(fi/M) acting on the particle, 

D2 ^ T) -F^). (1.1) 
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Formally F 11 is given by P,H 



F" = -i(<T + (2hlf. a - h^\ u ) u a u? , (1.2) 
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where u M = dz^/dr and is the tail part of the linear perturbation caused by the particle, i.e., the contribution 
from a part of the retarded Green function that has support inside the past lightcone. Thus our task is to evaluate 
the tail part of the Green function^. 

However, since the tail part depends non-locally on the background geometry, it seems practically impossible to 
evaluate it directly. As a way to avoid this difficulty, Mino proposed a very interesting approach || . He noticed that 
the direct part of the Green function, which has support only on the lightcone, is determined by the local geometry. 
Hence the corresponding metric perturbation can be evaluated solely with knowledge of the local geometry. On 
the other hand, in the black hole background, the linear perturbation is known to be described by the Teukolsky 
formalism ||, and there exists a method to obtain the general solution in a systematic manner |?]]. Then, using the 
transformation obtained by Chrzanowski ||, one can construct the corresponding metric perturbation. Hence, on the 
one hand, one may construct the full metric perturbation by solving the Teukolsky equation, and on the other hand, 
one can calculate the direct part of the metric perturbation by analyzing the local geometry. One can then obtain 
the tail part by subtracting the direct part from the full metric perturbation. 

This approach was further advanced by Mino and Nakano |J, and they succeeded in obtaining some interesting 
results. Yet, their method seems rather cumbersome. It involves many steps of lengthy calculations, mainly resulting 
from the different choice of gauges in which the full metric perturbation and its direct part are evaluated. 

In this paper, we propose a new, possibly more practical method to evaluate the tail part of the perturbation. 
As noted by Mino H, what one needs to do is to extract the direct part of the Green function. But this does not 
have to be done on the level of the metric perturbation. Instead, if one considers the Teukolsky equation in the 
spacetime domain, one can evaluate the direct part of the Green function for the Teukolsky equation directly, without 
transforming it to the metric perturbation. By Fourier-harmonic expanding the direct part obtained in this manner, 
one can then systematically evaluate the tail part. Then the corresponding metric perturbation and the reaction force 
can be calculated by using the Chrzanowski transformation M . 

A technical problem in this approach is that the Teukolsky equation has a very non-trivial angular spin weight 
dependence, even in the limit of flat spacetime, M — > 0. Thus the evaluation of the direct part becomes very 
involved. As a way to avoid this problem, we consider the transformation of the Teukolsky equation with spin weight 
s to a Klcin-Gordon-type equation (i.e., with spin weight s = 0). In the case of the Schwarzschild background, 
this is essentially the Chandrasekhar transformation that transforms the Teukolsky equation into the Regge- Wheeler 
equation In the case of the Kerr background, there is no known transformation of this kind. One possibility 

is to use the transformation given by Sasaki and Nakamura [fTl| in a perturbative manner, with respect to the spin 
parameter a/M of the Kerr geometry. In this paper, however, we focus on the case of the Schwarzschild black hole 
background and defer the extension to the Kerr case to a future work. 

There are also a couple of issues to be resolved before applying the present method to the actual calculation of 
the reaction force. First, the Teukolsky equation, by its nature, does not give the £ = and 1 spherical harmonic 
components of the perturbation. Hence the tail part of these components must be calculated separately. Second, and 
most importantly, since the reaction force is a gauge-dependent notion, its physical effect can be clarified only 
after we obtain the O(fi) correction to the gravitational radiation emitted to infinity, but what we have at hand is 
only the linear-order metric perturbation. Thus, strictly speaking, there is no point in evaluating the reaction force 
unless we have a method to evaluate the corresponding correction to the gravitational radiation. 

There have been other proposals for the regularization of the self force. Barack and Ori |l2| proposed a mode-sum 
regularization scheme. The calculation of the reaction force in their approach has two steps. First, each Fourier- 
harmonic mode of the bare self force is calculated. Then the sum over all modes is made subject to certain regu- 
larization procedure, which requires several regularization parameters. The calculation of the bare force on a scalar 
charge is done in |l3|-|T5|] . The analytical determination of the regularization parameters is carried out using a local 
perturbative analysis in ]l^Jl6[| . A method similar to that in [[l2| was proposed independently by Lousto ]l7j| . Very 
recently, using the mode-sum regularization prescription, the regularization of the scalar and electromagnetic self force 



in the spacetime of spherical shells has been discussed by Burko, Liu and Soen |18 |. 

This paper is organized as follows. In Sec. [b], we consider the Chandrasekhar transformation in the spacetime 
domain and discuss the nature of the direct part and the tail part of the Regge- Wheeler Green function. Then we give 



the Fourier-harmonic expanded form of the Teukolsky and Regge- Wheeler Green functions. In Sec. |I|, we formulate 
a method to calculate the tail part of the Regge- Wheeler Green function. Then, the tail part correct to the first 
post-Minkowskian (1PM) order is given under the short distance approximation, and an explicit expression of the tail 



in the first post-Newtonian (1PN) order is given. Then in Sec. IV, the procedure to obtain the regularized reaction 
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force is discussed. Finally, in Sec. [V| we summarize the results and discuss a possible way to resolve the remaining 
issues mentioned above. 



II. REGGE- WHEELER AND TEUKOLSKY GREEN FUNCTIONS 



A. Chandrasekhar transformation in spacetime domain 



The Teukolsky equation in the spacetime domain is expressed as 



-7 — ) o t + (r - M)d r + i—^-o v 
/ sm 



1 — s cot 



(2.1) 



where s = 0, ±1 or ±2, and Ak = r(r — 2M). We note that s \& has spin weight s with respect to the spatial rotation 
of the tetrad, but it is a spacetime scalar. The quantity S T is the source term whose explicit form can be found in ||. 
Although our interest is in the gravitational perturbation, and it is described by the Teukolsky equation with \s\ = 2, 
we leave the value of s unfixed until we need to fix it. As we see below, this makes it easy to identify the tail part of 
the Green function. 

The difficulty of directly dealing with the Teukolsky equation comes from the non-trivial angular dependence in its 
zeroth and first derivative terms. In particular, this causes non-conventional behavior of the Green function even in 
the flat spacetime limit. To overcome this difficulty, we consider a transformation of the Teukolsky equation. Since 
there is a one-to-one correspondence between the Teukolsky equation with spin weights s and —s (the Starobinsky- 
Teukolsky identities), we concentrate on negative spin weights s = 0, — 1 and —2. For s < 0, we introduce a new field 
variable S X by setting 



where 



s* = 9 



Ak 

r 



'dy 



-dt 



de + —. — -d v — s cot ( 



sm( 

i 

sin( 



■cL + s cot ( 



(2.2) 



(2.3) 



It should be noted that 9(9) operating on a quantity of spin weight s yields a quantity of spin weight s + 1 (s — 1) 
In particular, we have 



9 ( s Y em ) = [£(£ + 1) - s(s + 1)] 1/2 s+1 Y lm , 
9 ( s Y tm ) = - [1(1 + 1) - s(s - 1)] 1/2 s ^Y ern . 
99 ( s Y tm ) = - [1(1 + 1) - s(s + 1)] s Y im , 



(2.4) 



where the s Yg m are the normalized spherical harmonics of spin weight s. This implies that S X has spin weight 
s + \s\ = for s < 0. If expanded in terms of e~ lut and s Ye m , the radial part of the above transformation is just 
the Chandrasekhar transformation that relates homogeneous Teukolsky functions to Regge- Wheeler functions, apart 
from the overall normalization constant that depends on £ and to. A conceptually important difference is that here 
we do not associate S X with the spin weight s, but with s = 0. 



Inserting ( |2.2| ) into (2.1), we find 



sBs^ = 9' 



Ak 
r 



Or 



2M 



A K 



d t 



• 1*1+2 



2M 



-47r.,r. 



(2.5) 



We see that the equation for S X has the standard Klein-Gordon type structure with a potential term that vanishes 
sufficiently rapidly at infinity irrespective of the spin weight. Note that its radial part is just the Regge- Wheeler 
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equation. It then follows that S X can be naturally expanded in terms of the scalar (i.e., s = 0) spherical harmonics 
even for nonzero s, in accordance with th e di scussion above. A very important point to be noted is the following. If 
one constructs the Green function of Eq. ( |2.5| ) for S X, its direct part is exactly the same for all the spin weights. This 
is because the local causal structure of a hyperbolic equation depends only on the derivative terms of the equation 
(i.e., the V^V^ part). This implies that the part that should be discarded from the Green function is independent of 
s if done for the variable S X. In other words, any s-dependent terms in the Green function for S X belong to the tail 
part. 



B. Green functions for s ^ and S X 



We first consider the Green function s G rct for Eq. (2.5). It satisfies 



VV M + s 2 



2M 



5{t-t')S(r -r')S 2 (n-n r ) 



(2.6) 



Its Fourier-harmonic expansion is expressed as 



G rct (x,x') = I ,G& lw (p 1 r')oy< Bl («)ol r to.(fl')e- iw(, - t) ; 



s Gt mu> (r, r') 



1 X7(r>)Xf(r<) 



\u>\W{xy,X?) 



(2.7) 



where r> = max{r, r'}, r< = min{ r,r'}, oYi m is the usual (spin weight 0) spherical harmonic function, W(X^ p 1 X^ n ) 
is the Wronskian, 



_d_ 

dr* 



(2.8) 



and X^ p and X" 1 are the homogeneous solutions of the Regge- Wheeler equation having the asymptotic behaviors,^ 

' Ale' tur * + •Ag?e <wr * (r -» oo), 



e 



U lu> C 



(r -> 2M), 

(r -> oo), 
+ V ou±^r' ( r _>2Af). 



(2.9) 



It then follows that W = 2iw.A£, = 2iuV^. 



To 1PM order, the necessary formulas for the radial functions X 1 / 1 and X" p may be found in Poisson and Sasaki 
]20| and Leonard and Poisson plj . (Their explicit expressions are given in Appendix |A|.) For u> > 0, they take the 
form 



inz)^[j £ (z) + £ ^(z)+0(e 2 )], 
X7(z)= z \hf\z) + e$ p (z) + 0(e 2 ) 



(2.10) 
(2.11) 



where e = 2Mw, z = wr, and the functions Q n and £" p are defined in Eqs. (Al) and (A2). As noted in Appendix [A|, 
for oj < (i.e, for z < 0), they should be replaced by A™(|z|) and A" p (|z|), respectively. With this understanding, 
we find 



G imul {r,r') =«j j,(z<)^ 1) (z>)+e (j<(*<)£ p (z>) + CTi^K' ( Z> ) ) + 0(e 2 ) 



(2.12) 



Hence 



^In this paper, we follow the notation of Chrzanowski and Misner [n9| for the homogeneous solutions. 



4 



sGi niuJ Y£ rn (Q)Y£ 7n (Q / ) 

m 

(2£+ 1) [j,(z < )^ 1) (z > ) + e (jXz<)£ p (*>) +er(^<)^ 1) (^>)) +0(e 2 )] P/(/i) , (2.13) 



47T 



where /i = $7 • $7'. 

Now, given the Green function for S X, the Green function for is given as follows. The Teukolsky functions S -R™, 
s i?" p and the Regge- Wheeler functions S A™, S X^ P are related by the Chandrasekhar transformation, 

pin _ in f~i yin nup _ up ^ yup (9 

s-n,£ — sXl s^u s^-i i s-n,^ — sXf S^W S^-l J I* 2 - 1 v 

where S C W (s = 0, ±1,±2) are the Chandrasekhar operators |l^j2^,^l| and sxf 1 and s x" p are certain normalization 
coefficients (the case of s — is trivial, since qR = oX). For s = —2, the operator _2Cw is defined by 

C u =u»*f£f- 1 £r, (2.15) 

where 

d - _ / 2M\ 



C=f—+iu, f= 1 , (2.16) 

ar \ r / 

and we have omitted the spin index of -iC^ for simplicity and denoted it simply by C u , since we are interested only 
in the case s = —2. The coefficients ~2XT and -2X1 P are given by 

in = 16(l-2iA/o;)(l-4iMa;)(Ma;) 3 up _ 1 
Xf (*-l)*(*+l)(* + 2)-12»Afw ' Xl 4' 1 j 

where we have also omitted the spin indices for simplicity. In addition, in the rest of the paper, we denote by 
-04, in accordance with the conventional notation for the Weyl scalar. 

The Green function for the Teukolsky equation is constructed from the above homogeneous functions as given by 
Eq. (4.1) below. Except for the Wronskian W(R^ P , R™), however, the explicit form of the Green function for the 
Teukolsky equation is not needed for the following calculation, as will be clarified later. 



III. TAIL PART AT SHORT DISTANCES 



As usual, we write the Green function in the Hadamard form, 

s G rc \x, x') = 0(E(a), x') (^S(a) + 6(-<r) ,«) , (3.1) 

where a is half the squared geodetic distance between x and x' , u is a bi-scalar that depends only on the local 
geometry, and s v is a bi-scalar that describes the tail effects. We have a unique solution for each of these bi-scalars. 
As noted in the previous section, the direct par t of s G ret is independent of s, and the bi-scalar u is given simply by 
that in the case of the scalar d' Alembertian [E2| , 



det ler. 



a;0' 



y/g{x)g{x') 



1/2 

! (3.2) 



Let us denote the direct part of s G Ict by G d 



G d "( Xl x') = fl(E(s), x') 5(a). (3.3) 

47T 

Our task is to evaluate G dlr in the Fourier-harmonic expanded form at short distances. 

Since the full evaluation of the direct part is difficult, we consider its post-Minkowskian (PM) expansion. To do so, it 
is more convenient to work in the isotropic coordinates than in the Schwarzschild coordinates. The radial coordinates 
of these two coordinates are related as r = r/(l + M/2ri) 2 . Here r denotes the conventional Schwarzschild radial 
coordinate and rj denotes the isotropic radial coordinate. It should be noted that the isotropic coordinates satisfy 
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the harmonic gauge condition to 1PM order. The bi-scalar a(x, x') accurate to 1PM order can be found in a classic 
paper by Thorne and Kovacs [^3[ , or in Leonard and Poisson Q . 

We will come back to the explicit evaluation of the 1PM tail part later. Here, we first describe the formal procedure 
t o o btain the tail part. Let us express the Fourier- harmonic components of the full Regge- Wheeler Green function 
(0) as 

s G(£muj; x, x') = a G lmu {r, r') oY tm (n) Y tm (M) e -*"(t-0 . (3.4) 

Correspondingly, the Fourier-harmonic components of the direct part, which is s-independent, may be expressed as 

G diT (£moj; x,x') = Gf* w (r, r') Y im (Q) Y lm (n') e -*-(*-* / ) . (3.5) 

The tail part of the Green function, i.e., the regularized Green function, is obtained by subtracting the direct part 
from the full Green function: 

s G tail (emiu; x,x') = ,G^(r, r 1 ) Y tm {a) Y tm {Q,') e -M*-*') ; 

sG t hH UJ {r,r')= s Gtmw{r,r') — G^ u (r, r') . (3.6) 

Since this subtraction is done on a "mode-by-mode" basis, we call our procedure the "mode-by-mode regularization" 
of the Green function. 

Now let us return to the evaluation of the first post-Minkowskian (1PM) order tail. To 1PM order, the geodesic 
connecting x and x' can be approximated by a straight line on the flat background f23| , pT|| : 

^(A) = x' ft + XX> M ; X* = x tl -x'L l . (3.7) 

Here the affine parameter A is normalized so that £ M (0) = x lfJ - and £ M (1) = x^. We then have 



<j(x,x') = ~ (n^ + J h,j,v(€(X))d\ 



X^X 1 



-(-At 2 (l-F)+R 2 (l + F)) ; F = j ---dX, (3.8) 



2 v v ' v " ' 7o MX) 

where At = t — t' and R 2 = Sij(x 1 — x n )(x : ' — a;'- 7 ). The 1PM correction factor F can be calculated exactly as 

_ 2M (77 + r' I + R) _ 2M y, 2 ( R 2 " 

R (r T 4- r' — V}\ r, 4- r' ^ 



R 

AM 4M ( R 



(rj +r'j-R) n + r'j ^ 2n + 1 \n + r\ , 
2 

+ 0{R 4 ), (3.9) 



rj + r'j 3(r/ + r^) \ri+r'j 

where the short-distance approximation is applied to the second line. The bi-scalar u accurate to 1PM order is unity, 
because of the Ricci flatness of the Schwarzschild geometry. Thus the direct part of the Green function becomes 

G diT (x, x') = 25 {-At 2 (l — F) + R 2 (l + F)) 



4ir 

S(At-R(l + F)) 
l~^R 



(3.10) 



To subtract the direct part from the Green function, we take the Fourier-harmonic expansion of the direct part. 
The Fourier transformation of G dn is 



G A j{x,x') = J G d '"e iujAt dAt 



e iu)R{l + F) 



AnR 

Using the relation 

~~R~ 



AnR 

(l + iujRF + 0(M 2 )) . (3.11) 



Y,^^+^M^ri < )h ( p(iur I> )P e ( f i), (3.12) 



G 



and re-expressing rj in terms of the Schwarzschild radial coordinate r, we find 



Mz < )h ( l 1) (z > ) + 



— + — )M*<)f$\z>) 



+— hf ] (*>) + ^ n(z<) h^' £ ( Z> ) +0 (e\ \z - zf ) 



Pt(n) . 



(3.13) 



where z = u>r, z' = lot' and e = 2Mu, as before. Thus the 1PM direct part of the radial Green function is found as 



+ - fc) ^^>) + - 3t{z<) h (1 K(z>) +0 (e 2 , \z zf) 
z > z < / 



(3.14) 



The 1PM tail is obtained by subtracting the above from the full 1PM Green function, (2.12). 

Further performing the post-Newtonian expansion, the 1PN direct Green function is obtained explicitly as 



1 1 rr 



21+1 r' Vr'> 

M(£ + l) 1 uj 2 r' 2 M£ 1 u; 2 r 2 



x 1 + 



2 21- 1 r 2 2£ + 3 



+ 0(|r-r'| d ) 



(3.15) 



Expanding the full 1PM Green function (2.12) to 1PN order, and subtracting the above direct part (3.15) from the 
result, we find 



s 2 1 / r 



M 



M 



2i+lr' \r'J V {£+l)r' tr 



(3.16) 



It is worthwhile to note that the 1PN tail comes solely from the s-dependent part of the Regge- Wheeler equation. 
Thus the scalar d'Alenbertian V^V^ contains no 1PN tail, but only the correction to the light cone. In particular, 
this implies the 1PN tail is absent for a particle with scalar charge, in accordance with the recent result of Burko, Liu 
and Soen ||. 

IV. REGULARIZATION 

In this section, we give the procedure to calculate the regularized Weyl scalar, the tail part of the metric perturbation 
and the gravitational reaction force at the position of a point particle. 



A. Tail part of the metric perturbation 



From the homogeneous solutions X™ and X™ p of the Regge- Whe eler e quation, we obtain the corresponding solutions 
of the Teukolsky equation by the Chandrasekhar transformation (2.14). Using the method described in the previous 
section, the radial part of the regularized Green function for the Weyl scalar s Glmu>(T, r') is found as 



ifi/m»(r/) = -2 7r 



VF(i?; p , Rf) 



(4.1) 



where s = — 2 and C w (r) is the Chandrasekhar operator defined by Eq. (2.15). The factor — 2tt is inserted in the 
above formula for s Qimuj, so that the Weyl scalar -04 is given in the form 



2^ 



E 



dr'_ 2 ^ mw (r,r')T^ m (r')r'- 2 (r' _ 2 M) 



(4.2) 



with the source term T w f m defined in Appendix [5| in accordance with the definition of Poisson and Sasaki 
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Once we have the regularized Green function for the Weyl scalar, the regularized metric perturbation can be 
constructed by the Chrzanowski transformation ||. In the present case, we calculate the tail part (i.e., the regularized 
part) of the metric perturbation as follows. We decompose the radial part of the regularized Green function for the 
metric perturbation into parts having different spin weights as 

Gf^)(r,r') = ^gl sdu (r) sl Dl(r')^g lmuJ (ry), (4.3) 



where the coefficients s pi, and the operators s d w and s'D^ are defined in Eqs. (B13), (B14) and ( B15| ), respectively, 
of Appendix 0. Then the tail part of the metric perturbation is expressed as 



h^(£mu;;x) = [y v dr'r'- 2 (r' - 2M)- 2 G^(r, r') s T utm (r')^ Y em (9,, 

~{l^m v + rh^l v ) J 2M dr ' r '~ 2( - r ' ~ 2Af)- 2 G^)(^^) s ^ m (r')^ iY lm (9, <t>)e~ wt 
+fh^ j I2 [ dr ' r '~ 2{r ' ~ 2M y 2G \ls)^y)^^')] 2Y tm {9y)e- Mt ) : (4.4) 



where {l a } = { — 1, rj (r — 2M ), 0, 0} and {fh a } = {0, 0, 1, — i sin^jr /V2 are the Kinnersley null tetrad, and s T u h n is 
the complex conjugate of s Tui m - We note that the above metric perturbation is given under the ingoing radiation 
gauge condition, defined by 

V iM = °. V = °- ( 4 - 5 ) 

It should be noted that the Chandrasekhar transformation and the subsequent Chrzanowski transformation in 
their original forms can be applied only in the formal sense in the above procedure. In practice, one should modify 
those transformations as discussed in Appendix ^|, in order to subtract out the direct part induced by application 
of the original transformations. Let us explain the reason. The Chandrasekhar transformation involves the second 
order derivatives. If one takes the second-order derivatives of the Green function for the Regge- Wheeler equation 
in the spacetime domain, the step function 0(—a) in front of the tail part is transmuted to the direct part of the 
Green function for the Weyl scalar. This phenomenon shows up in the Fourier-harmonic domain, as demonstrated 
in Appendix Thus one should modify the Chandrasekhar transformation so as to subtract this induced direct 
part. As the Chandrasekhar transformation, the Chrzanowski transformation also involves the second-order radial 
derivatives. Hence, it is necessary to introduce a modified Chrzanowski transformation to subtract the induced direct 
part. In the actual calculation of the metric perturbation (or the reaction force), however, we do not introduce the 
modified Chandrasekhar and Chrzanowski transformations individually, but construct operators that contain only the 
first-order derivatives and that directly transform the Green function for the Regge- Wheeler variable to that for the 
components of the metric perturbation (or the reaction force) as illustrated at the end of Appendix ^|. The explicit 
expressions for the transformation operators are given in Appendix [b|. We note that thus constructed operators are 
same as the original operators with higher derivatives when applied to the full Green function. 

When we take the limit in which the position coincides with the particle, the following facts should be kept in 
mind. Let us denote the position of a particle by {z a } — {to, ro, 6*o, </>o} and a field point by {x^{ — {t,r,8,(p}. 
Because of the causal nature of the tail part, strictly speaking, the coincidence limit has to be taken by first taking 
the spatial limit as r — > ro, and then taking the limit t — > to- On the other hand, the coincidence limit for the angular 
coordinates may be taken at any time, as the metric perturbation depends on the angular coordinates only through 
the regular functions s Yi m - But as long as we employ the post-Newtonian (PN) expansion, we may take any order of 
the coincidence limit, because the light cone flattens out in the PN expansion. Therefore, when we evaluate a quantity 
in the coincidence limit, we take the limit t — > to first, followed by the spatial coincidence limit within the context of 
the PN expansion. After taking the limit, we perform the summation over £, m and lo. 

Here it is worthwhile to note the following. If we take the mode summation over i, m and u> before we take the 
spatial coincidence limit, the result depends on whether the radial limit is taken from r smaller than ro or larger than 
ro ■ This is a reflection of the fact that the first-order radial derivatives remaining in the final operators that transform 
the Regge- Wheeler variable to the metric perturbation (or to the reaction force) induce a 'quasi-direct' part that has 
a direction-dependent limit. Nevertheless, the avarage of the two different limits turns out to be the same as the limit 
obtained by first taking the coincidence limit before the mode summation. This can be explained by considering the 
first derivative of the step function in the spacetime domain. It has a direction-dependent limit, but this difference 
disappears when averaged over all directions. 
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B. Force for a circular orbit 



The form of the reaction force which describes the deviation from the geodesic motion on the background is given 
by Eq. (1.2). Here we focus on circular motion. Since only the t- and ^-components of the four velocity are non-zero 
in this case, the radial component of the reaction force can be expressed as 



F r {x) 



1 



(2h rt - t - h t t;r) "V 

- (2/l r ; — /l00 ;r ) U^U 4 " + 2 (htr;<t> + hffyrtf 



tail 



(4.6) 



where [• • -] tail indicates that h^ v:a should be interpreted as the tail part h^: a . Specializing to the case of Schwarzschild 
background, we have 



F r (x) 



2M/, , , , 

2h rt ,t - h t t,r -^—hrr ) U U 

+ {2hr<i>,<i> ~ h<p<p. 
+2 [htr,4> + h^ r ^t 



j- 2rsia 2 9fh rr ) u^u 4 



tail 



(4.7) 



Since there is no effect of the gravitational radiation at 1PN order, we may consider only the radial component of the 
reaction force to that order. Physically, the radial component describes the correction to the radius of the orbit that 
deviates from the geodesic in the unperturbed background. 

As discussed in the previous subsection, we c onst ruct the operators that directly trasform the Regge- Wheeler 
variable into the reaction force expressed in Eq. (4/7). At 1PN order, we find that only the h^ 1 ]. term contributes. 
Since this term involves the operation d r oA. u , where oA w is defined in Appendix [b|, Eq. ( B19Q , we reduce the second 
order derivative in it by using the homogeneous Regge- Wheeler equation, as discussed in the preveous subsection. 
The result is 



d r0 A u 



( \„ /r J 2M\ ( 9M\ 1,, M ( 3M 
= f --£(1+1) ( 1 ) -iwrf 1- — ) + -cjV + — ( 1 



4M 2 



r J 
1 + 2 i ujr - 

2 e 2 + 2 1 + 3 



2r J 2 



r J 
6M 
r 



+ 2w V 



+ 2 i lu \2t + 21 + 5 



r \ r 
12M\ 



~d r 

r 



r j 



(4.8) 



This operator is used to derive the 1PN reaction force term. 

After summing over £, m and uj, and taking the coincidence limit, the gravitational reaction force to 1PN order is 
found to be 



F r {z) 



11 Mn 
3rg 



(4.9) 



We warn that this result by itself has no physical significance. First, it contains contributions only from harmonics 
with £ > 2. We have to calculate the £ = and 1 components to make the force complete. Second, it is obtained in 
the ingoing radiation gauge, and since we lack a method to calculate the gravitational radiation to second order in /i, 
the only thing we can do at present is to compare the result with the standard 1PN formula given in the literature. 
To do so, we have to make a gauge transformation from the ingoing radiation gauge to the harmonic gauge. Only 
after this procedure can we determine if our result is valid. 



V. DISCUSSION 



In this paper we have proposed a new method to derive the gravitational reaction force on a point particle in 
the Schwarzschild background. In this method, the tail part of the metric perturbation, which is responsible for 
the reaction force, is obtained by first considering the Fourier-harmonic expansion of the Green function for the 
Regge- Wheeler equation. After calculating the tail part of the Green function, the result is transformed to the 
metric perturbation through the Chandrasekhar transformation followed by the Chrzanowski transformation. Since 







the extraction of the tail part is done at the level of the Green function in the Fourier-harmonic expanded form, 
we call our method the mode-by-mode regularization of the Green function. An advantage of our method is that 
it is relatively easy to obtain the tail part of the Green function for the Regge- Wheeler equation if we adopt the 
post-Minkowskian (or post-Newtonian) expansion. 

However, there are still a couple of issues to be resolved. One is that the £ = and 1 components are not included 
in the metric perturbation obtained with the Chrzanowski transformation. Hence the tail part in them must be 
derived by some other method. The £ = and 1 components of the full metric perturbation can be solved exactly 
in the Regge- Wheeler- Zerilli gauge, as shown by Zerilli ]2^| . Therefore, what we need is a method to identify the 
tail part in the £ = and 1 components. The other, most crucial issue is that the physical meaning of the resulting 
reaction force remains unclear, since the reaction force is a gauge-dependent notion. For an ultimate resolution of 
this essential problem, we need to develop a second-order perturbation theory so as to make it possible to calculate 
the 0(/i) correction to the gravitational radiation. 

Nevertheless, at low post-Newtonian orders, there is a way to resolve these issues. The equations of motion for 
a two-body system and the gravitational radiation from such a system are investigated extensively in ]25|-p8| in the 
harmonic gauge. Hence, by performing the gauge transformation to the harmonic gauge, we should be able to clarify 
the physical effect of our result, as well as to examine the consistency of our result with those of these previous works. 
As a first step, we are currently working on a consistency check at the 1PN order |29[| . 
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APPENDIX A: THE HOMOGENEOUS 1PM REGGE- WHEELER FUNCTIONS 



Here we recapitulate the formulas for the homogeneous 1PM Regge- Wheeler radial functions X 1 / 1 and A" p , which 
can be found in Poisson and Sasaki |2(J and Leonard and Poisson [gTJ . 
The formulas are 



Xf{z)= 2A i £ J (- l ) e+1 e^ 2 ^\l + | e) z [j t (z) + e&(z) + 0(e 2 )] . 



&(z) = -A e (z)j e (z) + B e (z)n e (z) - 



2£(2£+l) 



2(£ + 1)(2£ + iy 



A7 P (z)= i<+i e -fc(M2 e )-/W (1 _| e) , ^(D (z ) + e ^p + (e 2 ) 



where e = 2Mlu, z = tor, and the functions A^{z), Bi(z) and Ct(z) are given by 



Ai{z) = Si(2z) + z 2 n { 



z 2 n p {z)j p {z) , 



B e (z) = Ci(2z) - 7 - ln(2z) + z 2 j (z) 2 + ^ (- + J-) z 2 Jp (z) 2 , 

C t (z) = 2 (| - Mz) + i (B t {z) + 7 + ln(2z))) - i (l + £ (± + -L^j z{-iy 



{*)), 



(Al) 



(A2) 



(A3) 



for i > 1, and 
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A (V) = Si(2z), 
B Q (z) = Ci(2z)- 7 



ln(2z) 



Co (2) = 2 ( - - Si(2«) + / 



iCi(2z)) 



(A4) 



for £ — 0. It should be noted that the expressions (Al) and ( |A2| ) are valid only for w > 0. For u < (i.e, for z < 0), 
they should be replaced by X™(|z|) and X^ p (|z|), respectively. 



APPENDIX B: THE CHRZANOWSKI TRANSFORMATION 

In this appendix, we summarize the Chrzanowski method Q to construct the metric pertubation from the Weyl 
scalar and describe the modification of differential operators necessary to extract the tail part of the metric pertur- 
bation. 

Using the method of Chrzanowski [^| , we can generate metric perturbations in either the in-going radiation gauge 
or the out-going radiation gauge from the Weyl scalar. We first introduce the following differential operators on the 
Kerr spacetime: 

h™ n = -(5 + a + 30 - f)(5 + + 3t) , 
K m = -(D-p + M-e)(D + 3p + 4e) , 
t£ m = -1/2{(D + p-p + e + 3e)(S + 4/3 + 3f) 

+ (5 + 3p-a-TT-f)(D + 3p + 4e)} , (Bl) 



and 



hff* = -(<$ + 7f - 3a - (3){8 - 3tt - 4a)( j o)" 4 , 

\-4 



hT m = -(A + 7 - 37 + p)(A - 47 - 3p)(py 
hfg = -l/2{(5 + r + 7f - 35 + /3)(A - 4 7 - 3/2) 

+ (A - 7 - 3 7 + p - p)(<5 - 3tt - 4a)}(p)~ 4 . (B2) 

Here Z? = Z M V M , A = ri^V^, <5 = m^V^ and 5 = m^V^, where n^ 1 , m/*, m* 1 } is the Kinnersley null tetrad, p, /i, 
r, 7r, e, 7, a and (3 are the spin coefficients of the Newman-Penrose formalism |30f] , and the overline indicates complex 
conjugation. 

With these differential operators, we have the following metric perturbation in the Fourier-harmonic expansion in 
the ingoing and outgoing radiation gauges: 

h%{£mw) = {l^luKn + fh^hZm - (i P *n„ + rn^)h™ m }„ 2 R emL ,(r) 2 Y ern {9, fie** , (B3) 
<(H = {n^ftjf + m^fc^ - (n M m„ + ™ A1 n 1 ,)^ t } 2 i?, mLJ (r)_ 2 r, m (0, 0)e-*"* . (B4) 

Here s i?; mw (r) is the Teukolsky function of spin index s, and s Ye rn (9, <fi) is the spin weighted spherical harmonic 
function. 

Now we specify the background spacetime to be a Schwarzschild blackhole. We assume operands are expanded in 
the Fourier-harmonic functions. We use the pure radial and angular operators introduced by Chandrasekhar [fTo) , 

„ „ ir 2 bj r — M „ ir 2 oj r — M . . 

V n =d r - — + 2n— , 2?t = d r + —— + 2n— , (B5) 

£■<% = dg + m + n cot 9 , C' n = d$ — m + n cot 9 , (B6) 

sm 9 ' sin 9 

where Ak = r[r — 2M). The Newman-Penrose operators and spin coefficients become 



D = V , A = --4x>, f 



2r 2 »' V2 



T/l CB7) 



1 cot0 A K M . . 

' = - a = /3 = ^' ^ = -2^' 7= 2^' (B8) 
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The radial differential operator £ introduced in Eq. ( 2.16| ) is related to the operator D in the above as £ — fD, where 
/ = (r — 2M)/r. Then all the operators defined in Eqs. (Bl) and ( [B2| ) are expressed as 



tin _ _}_f f /,in 



1 ( m 2 \ 



.out 

n ll 



r 

——nn h ont 

2 i "'l i -'2 ' 'Hfh 



£2 



2?o 



2?o 



2\/2 V 



A 2 / 1 
£ t tout _ _'±K_( T) \ 1 
■°2 ' "'mm — ^ \^2 



(B9) 
(BIO) 



In the remainder of this appendix, operators used to derive the metric perturbation from the Regge- Wheeler variable 
are given. We consider the metric perturbation in the ingoing radiation gauge. 

First, note that, when oper atin g on a quantity of spin weight s, the angular operators £ s and £_ s are just the 
"cdth" operators given in Eq. (2.3), except for the over-all sign: 



9 = -£L. , & = -£.. 



(Bll) 



Hence 



Accordingly, we put 



£2 ( 2 Y em ) = ( 2 Y em ) = [(£ -1)(£ + 2)\ 1 ' 2 1 Y lm , 
C1JC2 ( 2 ^m) = & 2 ( 2 *im) = [(£ - !)£{£ + l)(£ + 2)] 1 / 2 r £m . 



C = oPi od-u , 

Kim = -lPi l<h 
h-mm = -2Vl 2 d u 



(B12) 



where the coefficients s pi (s = 0, —1, —2) are defined by 

o Pi = 2[(£-l)£(£+l)(£ + 2)} 1 / 2 , 
_ 1 ^ = 2[2(£-l)(£ + 2)] 1 / 2 , 

-2P( = 1 , 

and (s = 0, 1, 2) are the radial differential operators defined by 



(B13) 



1 



iL = --rf-'Cr- 2 , 

2L = -r- x r x £^r x Lr- 



(B14) 



Next, we rewrite the Green function for the Weyl scalar ipi in a more convenient manner, so that the source term is 
divided into the tetrad components of the energy-momentum tensor that have specific spin weights, s T^^ m , as defined 
by Eq. (D2). For this purpose, we introduce the operators S D W and S D£, defined in Poisson and Sasaki | [20[ : 

D U = Dl = r 4 , 
- 1 D u =r 2 f£r s f- 1 , 
-2D LU =rf£r 4 r 1 £r, 
- 1 Dt = -r 7 £r- 2 , 



-2D 



t - J> 



f£r i r 1 £r- 



(B15) 



Using these operators, we can rewrite the expression for ipi given by Eq. (4.2) as 
dui 



r ip4:(x) 



2tt 



E 

tin 



dr' - 2 g tmu (r, r') T ulm {r') r'- 2 (r' - 2M) 



2AI 



-2Ylm{tye 



dtu x -> 
2^^ 



J2 / dr'r'~ 2 (r' - 2M)~ 2 (s) G tmu} {r, r') s T^ ro (r') 



2A1 



-2Yt m (Q)e~ 



(B16) 
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where the Green Functions i s )Gi mul {r, r') (s = 0,-1, —2) are defined by 

(s)Gi mu (r, r') = s pi s Dl(r') -2Gimw(r, r') 



(B17) 



with s Dlj(r') acting on the r'-dependent part of -iQimuir, r') . The relation between the Teukolsky radial function 
-2-R{mu('') to be used in the Chrzanowski transformation (p33|) and the radial part of r 4 ip4(x) is 



32 



(oPi) 2 



V / dr'r'- 2 (r' -2M)- 2 {s) G tmuJ {ry) s T ujlm {r') 

„ JIM 



(BIS 



Here, however, we wish to construct a Green function that directly gives the metric perturbation. For this pur pose , 
we introduce the differential operators S A W = s d<_j C w and S A W = gD^ C w , where s d u and S Z?J, are defined by Eqs. (B14) 
and (B15), respectively, and C w is the Chandrasekhar operator defined by Eq. (2.15). Explicitly, S A W and S A W are 
given by 



iA„ = -\rf- 1 CfCf- 1 Cr = ^r" 1 /" 1 -i^, 
where the operators s r w are introduced in |[20[ and are given by 



-2J- u ■ 



(B19) 

(B20) 
(B21) 



o 



3M 



2 1 +»«r r/— + / \£{l + l) 



dr 



6M 



3M 

2iujr I 1 Hwr 

r 



-ir„ 



-f 



(£(£ + l) + 2iu;r)rf— + 1(1 + l)(f + i ur) - 2{ur) 2 
dr 



l T UJ = fCr i f- l Lr- l fCf- 1 Cr 



2 (I -!)(£ + 2) + 



6M 



r /4- + (t ~ W + 2)(£{i + 1) + 2iujr) - 12 f — 
dr r 



(B22) 
(B23) 
(B24) 



With these newly introduced operators, the regularized Green functions for the components of the metric perturbation 
can be written as 



Gffi)ir, r') = --^ ^V^) XT X7 M s >Pi s^{r) S -A w (r') _ 3 Ggl(r, r>) , 



(B25) 



where the coefficeints \™ and x" p are defined in Eq. ( [2.17 ). Note that the operators S A W and 8 T U defined above 
become the modified differential operators by reducing the higher r-derivatives using the homogeneous Regge- Wheeler 
equation. As discussed in Sec. IV, these operators transform the tail of the Regge- Wheeler Green function into the 
tail of the Green functions for the metric components. If the original higher differential operators are used, the direct 
part is partially induced, as shown in Appendix |c[ 



APPENDIX C: THE INDUCED DIRECT PART OF THE GREEN FUNCTION 

In Appendix |b], we introduced modified differential operators that differ from the original Chandrasekhar and 
Chrzanowski operators in order to obtain the tail part of the metric perturbation. This is because the direct part is 
partially induced by taking derivatives of the tail part of the Regge- Wheeler Green function, and such an induced part 
must be subtracted out. In this appendix, we explicitly show this phenomenon at 1PN order. Then we illustratively 
summarize our prescription to subtract the induced direct part from the Green function for the metric perturbation. 

The homogeneous radial Regge- Wheeler function satisfies 



0i 



Xt = o, 



(CI) 
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where z = wr, e = 2uM, and the operators Oq and 0\ are given by 



dz 2 



£{£ + !) 



O x -- 



2tf_ 

z dz 2 



1 d 

z 2 dz 



£{£ + !) 



(C2) 



The Green function is constructed from the homogeneous functions as X 1 g n (z)X £ p (z'). Therefore we consider this 
form. The homogeneous Regge- Wheeler functions can be expanded to 1PM order as 

X t = X + e (Xf + s 2 X{) + 0(e 2 ) , (C3) 

where Xf and X\ are the parts of the 1PM correction independent of s and dependent on s, respectively. When 
we further adopt the post-Newtonian expansion, they correspond to the parts of the homogeneous function that 
contribute to the direct and tail parts of the Green function, respectively, as discussed in Sec. [II. Specifically, the 
direct part of the Green function consists of Xq(z)Xq{z') and Xo(z)Xf(z') + Xf(z)Xo(z'), and the tail part consists 
ofX Q (z)X t 1 (z')+X t 1 (z)X (z'). 

We substitute Eq. (C3) into Eq. (CI) and collect terms of the same order in e. Then we find 



O X = , 

O.Xo + OoX^ + s 2 \O X\ + ± 

Since Xq and Xf are independent of s, it follows that 

O x X + O X* = 0, 

o x[ + \x = o. 



X = . 



(C4) 



(C5) 



We find from the second equation of the above that the tail part of the Regge- Wheeler Green function is transformed 
into the direct part by the operation Oq. 

To subtract out such an induced direct part, it is therefore nessesary to reduce the second-order derivative by 
the equation OqX{ = whenever second- or higher-order derivatives operate on the Regge- Wheeler Green function. 
Although we do not have a general proof, we expect the same to be true at any order of the post-Newtonian expansion. 
Thus, in general, the prescription is to reduce the higher-order derivatives on the tail part by the equation 



Al = 0, 



(C6) 



where X\ is the part of the homogeneous Regge- Wheeler funtion that constitutes the tail part of the Green function. 

It may be useful to summarize the above prescription illustratively. First, one constructs the full Regge- Wheeler 
Green function, say Grw- Then the full Green function for the metric perturbation, say G, is obtained as 



G(x,x') = 0(x)0'(x')G RW (x,x'), 



(C7) 



where O and O' are some higher-order differential operators. As long as x ^ x', one can then reduce the higher-order 
derivatives in O and O' by using the homogeneous Regge- Wheeler equation without any problem. Let us denote the 
reduced operators thus obtained by O and O' . Then we have 



G(x,x') = &(x)6'(x')G KW (x,x'). 



(C8) 



By construction, the operators O and O' contain only the first derivatives. Now we divide Grw into the direct part 
G^Jy and the tail part G^. Then the tail part of the Green function for the metric is given by 



G tail (x, x') = 0(x)d'(x')G^(x, x'). 
This is our prescription to obtain the tail part of the Green function for the metric perturbation. 



(C9) 



14 



APPENDIX D: NOTATION FOR THE SOURCE TERM 



The source term of the inhomogeneous Teukolsky equation is constructed from the energy-momentum tensor of a 
point particle: 

T af3 {x) = n [ dTu a u p 5 {i) [x - x'(t)\ . (Dl) 



The first step is to obtain the projections by the null tetrad, T a p n a n* 3 , T a p n a mP and T a p m a mP . Then we calculate 
the Fourier-harmonic components according to 

s T utm = dtdn s T{x) s Y lm {n)e 1 ^ , (D2) 

where the quantities S T (s = 0, — 1, —2) arc defined by 

T = T a0 n a n & , _ X T = T a(3 n a m 13 , _ 2 T = T a0 m a m? . (D3) 

Then the source term T^g m for the radial Teukolsky equation is obtained by operating on s Tuim with relevant 
diferential operators and then by summing over s as 

T u £ m = ^ s pe sD^ sT^im , (D4) 



where the a D u are the differential operators defined in Eq. (jBlq). 

We note that the source of the Teukolsky equation ( |2.1| ) is related to the above T w ^ m as 

-2%, im = ^- / dtdn_ 2 T{x) _ 2 y /m (n) e Mt = ^ . (ds) 

2ir J 2r z 

In this paper, we use the notation of Poisson and Sasaki [ p0[ for the source term and use T W £ m . 

Note that Chrzanowski || uses different notation for the source term T in the Teukolsky equation. Denoting 
Chrzanowski's T by Tc, it is related to T u i m and -2%j£ m as 

T c = -2 7T r w£m = -4 7r r 2 -2%im ■ (D6) 
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